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Abstract 

Unlike the classical polynomial case there has not been invented 
up to very recently a tool similar to the Bernstein-Bezier representa- 
tion which would allow us to control the behavior of the exponential 
polynomials. The exponential analog to the classical Bernstein poly- 
nomials has been introduced in the recent authors' paper [2], and this 
analog retains all basic propeties of the classical Bernstein polynomi- 
als. The main purpose of the present paper is to contribute in this 
direction, by proving some important properties of the Bernstein ex- 
ponential operator which has been introduced in [2|. We also fix our 
attention upon some special type of exponential polynomials which 
are particularly important for the further development of theory of 
representation of Multivariate data. 
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1 Introduction 



With the great increase in the number of businesses making a presence on 
the Internet and the increase in the number of cyber-customers, the chances 
of computer security attacks increase daily. How much should an enterprise 
or organization budget to defend against Internet attacks? 

As it is well know the processes modelled by ordinary differential equa- 
tions and the stochastic differential equations describe many important prac- 
tical events, and thus provide very efficient models for processes coming from 
physics, finance, etc. As is well known their solutions may be approximated 
efficiently by means of finite linear combinations of exponential functions, 
called sometimes exponential polynomials. Let us mention some recent ap- 
plications to exponential smoothing used for time series models which fit best 
for forecasting Internet Security attacks, [16] . 

Going further in this direction, if one considers a process which expe- 
riences jumps or unsmoothness, then it may be approximated by means of 
piecewise exponential polynomials. 

For the control and design of polynomials the famous Bernstein polyno- 
mials provide a very efficient way to control their behaviour. Their further 
development as the Bezier curves provides a very efficient method to design 
a prescribed form, which is an indispensable tool for Computer Design. In a 
similar way the theory of B— splines provides us with an indispensable tool 
for the control and design of the classical splines. 

Let us add to the above that for the purposes of fast recognition, repre- 
sentation and compression of curves and surfaces such tools as the Bezier 
curves and B— splines curves are very important. Since many of the real pro- 
cesses produce observational and surveillance data which carry exponential 
character, it is clear that it is very important to have an efficient tools for 
their representation, design and control. 

It is curious to mention that unlike the classical polynomial case there 
has not been invented up to very recently a tool similar to the Bernstein- 
Bezier representation which would allow us to control the behavior of the 
exponential polynomials. The exponential analog to the classical Bernstein 
polynomials has been introduced in the recent authors' paper [2J, and this 
analog retains all basic propeties of the classical Bernstein polynomials. 

The main purpose of the present paper is to contribute in this direction, 
by proving some important properties of the Bernstein exponential operator 
which has been introduced in [2]. We also fix our attention upon some 



2 



special type of exponential polynomials which are particularly important for 
the further development of theory of representation of Multivariate data. 

Let us recall at first shortly the concept of Bezier curves and its relation- 
ship to the Bernstein polynomials: Let bo,...,b n be vectors either in M 2 or 
M. 3 and t G R, and interpret b& as a constant curve, i.e. that b° k (t) := bt 
for k — 0, n. Then define new curves b\ (t) for k — 0, n — 1 by b\ (t) = 
{1 — t)b° k (i) + tb k+1 (t) . Repeating this process one arrives at new curves 

¥ k (t) = (1 - t) 6JT 1 it) + tbl~\ (t) for * = 0, n - r 

In the last step, i.e. for r = n, one finally obtains exactly one curve 

&s (*) = (i - 1) b-- 1 (t) + wr 1 ^ , 

the so-called Bezier curve. The polygon formed by b ,...,b n is called the 
Bezier polygon or control polygon. The Bezier curve has the property that the 
curve &o (t) is in the convex hull generated by the points bo, b n . Moreover 
the two points bo and b n are fixed, i.e. 6q (0) = bo and b$ (1) = b n . An explicit 
form for the Bezier curve is 

n 
A;=0 

where p n> fc (t) := (™)t fc (1 — are called the Bernstein basis polynomials. 
In the sequel we shall focus on generalizations of Bernstein basis polynomials 
which have arisen recently in Computer Aided Geometric Design for modeling 
parametric curves. Instead of the basic polynomials 1,2, ....,x n one consider 
different systems of basic functions f , f n , e.g. 

l,x,...,x , cosx, sins, 

which are better adapted to curves in spherical coordinates, see e.g. [B], [24|, 
[33] and [7] . In mathematical terms it will be required that the linear span of 
the basis functions f , f n forms a an extended Chebyshev system. Recall 
that a subspace U n of C n (I) (the space of n-times continuously differentiable 
complex-valued functions on a interval I) is called an extended Chebyshev 
system for a subset A C /, if U n has dimension n + 1 and each non-zero 
function / G £/„, vanishes at most n times on the subset A (counted with 
multiplicities). A system p n ^ G U n , k — 0, n, is a Bernstein-like basis for 
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U n relative to a, b G J, if for each k = 0, n the function p n> fc has a zero of 
order k at a, and a zero of order n — k at b. 

In the following we shall consider Bernstein basis polynomials for the 
space of exponential polynomials •£ , (a ,...,A») (induced by a linear differential 
operator L) defined by 



where Aq, A n are complex numbers, and L is the linear differential operator 
with constant coefficients defined by 



Exponential polynomials are sometimes called L -polynomials, and they pro- 
vide natural generalization of classical, trigonometric, and hyperbolic poly- 
nomials (see [31]), and the so-called P-polynomials considered in [28J. Expo- 
nential polynomials arise naturally in the context of a class of multivariate 
splines, the so-called polysplines, see [17], [20]). Important in this context 
are exponential splines associated with linear differential operators L 2s +2 of 
order 2s + 2 of the form 



which are parametrized by real numbers A and /x. According to the above 
notation A n = (Aq, A n ) we set n = 2s + 1 and define the vector 



containing s + 1 times A and s + 1 times ji. 

Let us return to the general theory of exponential polynomials, and let us 
recall the general fact (cf. [27J) that for A n = (A , A n ) G C n+1 there exists 
a unique function $a„ G £?(a ,...,a„) such that $a„ (0) = .... = ^ (0) = 
and $ A ™' ) (0) = 1. We shall call $a„ the fundamental function in En ,...,x n )- 
An explicit formula for $a u is 



E iXo ,..., Xn) := {/ G C°° (R) : L/ = 0} 



(1) 



L:=L(A) := f-f-Ao) .... 

\ax J \dx J 



(2) 





(3) 



(x) := [A ,...,A n ] e : 



(4) 



2m J Tr (z - A ) ... (z - A n ) 



4 



where [A ,...,A n ] denotes the divided difference, and T r is the path in the 
complex plane defined by T r (t) = re lt , t G [0, 2tt], surrounding all the scalars 
Ao, A n . The fundamental function $a„ is an important tool in the spline 
theory based on exponential polynomials (see [27] ) and in the wavelet analysis 
of exponential polynomials, see [8], [23], [T7] , [18], [T9] . 

In Section [2] we provide some basic results about the fundamental func- 
tion, in particular we derive the Taylor expansion of the fundamental function 
$(A ,...,a„)- Section [3] is devoted to the analysis of the fundamental function 
with respect to the eigenvalues ([3]), and in order to have a short notation we 
set using (TJJ 

$2s+i {x) := $a 2s+1 (a iA1 ) for A = 1 and /i = -1. 

We shall give explicit formula for $2s+i {%) • So far it turns out that the 
following recursion formula is more important: 

* 2 - +3 {x) = ^T^iTTf 2 - 1 {x) ~ tri* 2 - +1 w ■ 

We shall prove the non-trivial fact that for each fixed x > the sequence 

$2, 0) 



$2s- (x) 



(5) 



for s — > oo; here $2s {x) is the fundamental function with respect to the 
vector A2 S consisting of s+1 eigenvalues 1 and s eigenvalues —1, and similarly 
$2s- {x) is the fundamental function with respect the vector A2 S consisting 
of s eigenvalues 1 and s + 1 eigenvalues —1. 

In Section H] we shall determine some generating functions: we prove that 

X 

*2s+i (x) ■ y s = y=^j sinh (xy/y + l) 



s=0 



and that 



oo ^ 

$ s (x) ■ y s = —^^=^= sinh (x^/y 2 + 1 ] + cosh [x\/y 2 + 1 J . 

,=o VV + 1 v y v y 

In Section [5] we give a detailed introduction to the notion of a Bernstein 
basis p n ,k(x),k = 0, ...,n, in the setting of exponential polynomials with 
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arbitrary (complex) eigenvalues. We shall give a new proof of the result that 
there exists a Bernstein basis p ni k(x),k = 0, ...,n, in -E(a ,...,a„) f° r points 
a 7^ b if and only if -E(A ,...,A n ) is a Chebyshev space with respect to a, b. In 
the case that Ao, A n are real it is well known that S(A ,...,A n ) * s a Chebyshev 
space with respect to the interval [a, b] and in this case p n ^ (x) may be chosen 
strictly positive for x in the open interval (a, b). 

In Section [H] we shall derive recursion formulas for the Bernstein basis 
Pn,k (x) , k = 0, ...,n for the special system of eigenvalues A 2s +i which have 
been used in Section |3j 

A remarkable result was proved recently in [5], |26j for certain classes 
of extended Chebyshev systems U n : Assume that the constant function 1 
is in U n ; clearly then there exist coefficients k = 0, ...,n, such that 1 = 
Ylk=o a kPn,k, since p n ,k,k = 0,...,n, is a basis. The normalization property 
proved in [5] and [26] says that the coefficients are positive. It seems 
that the paper [2] addresses for the first time the question whether one can 
construct a Bernstein- type operator based on a Bernstein basis p n ,k-, k = 
0, ...,n, in the context of exponential polynomials, i.e. operators of the form 

n 

B n (x) := [B n f] (x) :=^2 a kf Pn,k 0) (6) 

k=0 

where the coefficients «o, ...,a n should be positive and the knots to, ■■■,t n in 
the interval [a, b] . In [2] the following basic result was proven: 

Theorem 1 Assume that Ao,...,A n are real and Ao 7^ Ai. Then there exist 
unique points t < t\ < ... < t n in the interval [a, b] and unique positive 
coefficients a , ...,a n such that the operator B n : C [a, b] — > Erx ,...,\ n ) defined 
by (TJP has the following reproduction property 

[B n (e x °-)] (x) = e Xox and [B n (e Al ')] (x) = e^ x . (7) 

The positivity of the coefficients ao, oc n is related to the above-mentioned 
normalization property. Theorem [T] says that property ([7D can be used for 
defining knots t , t n and weights a , a n for an operator of the form (jB]). 
In the classical polynomial case this means that the Bernstein operator B n 
on [0, 1] has the property that B n (l) = 1 and B n (x) = x for the constant 
function 1 and the identity function. 

It follows from the above construction that the operator B n defined by ([6]) 
and satisfying (j7j) is a positive operator. Using a Korovkin-type theorem for 
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extended Chebyshev systems the following sufficient criterion for the uniform 
convergence of B n f to /, for each / e C [a, b] has been given in [2J. Here 
we use the more precise but lengthy notation P(x ,...,x n ),k instead of p n ,k for 
k = 0, n. 

Theorem 2 Let Ao,Ai,A2 be pairwise distinct real numbers and let A n = 
(Ao, Ai, A n ) G R n+1 with possibly variable Xj = Xj (n) for j = 3, ...,n. For 
each natural number n > 2, and each k <n, define the numbers a (n, k) and 
b {n, k) as follows: 

a(n,k) :=hm P^M,...,x n lk (x) 
x ~* b P(MA 2 ,....,\ n ),k (x) 

b(n,k) . = i im P^MM,-..-M)A^) (g) 
x ^ b P(\i,\ 2 ,...,X n ),k (X) 

Assume that for n — ► oo uniformly in k holds 

t k (n) - t fc _x (n) — > 0, (10) 

where tk = tk (n) are the points determined by Theorem^ and for n — > oo 
uniformly in k holds 

i2i^l^A 2 -v (11) 

Then the Bernstein operator B^x ,...,x„) defined in Theorem^ converges to the 
identity operator on C([a,b],C) with the uniform norm. 

Theorem [1] applied to the system A 2s +i considered in Section [3] shows that 
there exist unique points t < t\ < ... < t n in the interval [a, b] and unique 
positive coefficients a ,...,a n such that the operator B n : C[a,b] — > E\ 2s+1 
defined by (jSj) has the property 



Ba 2s+1 (e x ) = e x and B A2b+1 (e x ) = e 



X 



The contributions in this paper may serve to investigate the question whether 
the Bernstein operator B\ 2g+1 converges to the identity operator for s — > oo. 
Since in this example one has only two different eigenvalues Theorem [2] has 
to be modified. Note that for k = n the coefficient a (n, k) defined in (jSJ) is 
equal to 

P(\ ,x 2 ,....,x n ),n (x) = <t> 2s (6 - a) 

x ^ b P(X 1 ,X 2 ,....,X n ),n (X) $2s- (b - a) 

and in (jSJ) we have proved that these numbers converge to 1. 



7 



2 Taylor expansion of the fundamental func- 
tion 



We say that the vector A n e C n+1 is equivalent to the vector e C n+1 if the 
corresponding differential operators are equal (so the spaces of all solutions 
are equal). This is the same to say that each A occurs in A n and A' n with the 
same multiplicity. Since the differential operator L defined in (j2J) does not 
depend on the order of differentiation, it is clear that each permutation of 
the vector A n is equivalent to A n . Hence the space E(x ,...,\ n ) does not depend 
on the order of the eigenvalues A , A n . 

We say that the space Erx ,-,X n ) is closed under complex conjugation, if 
for / e -E(A ,...,A n ) the complex conjugate function / is again in Erx ,...,x n )- It is 
easy to see that for complex numbers A , A ra the space Er\ \\ is closed 
under complex conjugation if and only if there exists a permutation a of the 
indices {0, n} such that Xj = A CT (j) for j = 0, n. In other words, £7(a ,...,a„) 
is closed under complex conjugation if and only if the vector A n = (Ao, A n ) 
is equivalent to the conjugate vector A n . 

In the case of pairwise different Xj, j = 0, ...,n, the space £?(a ,...,a„) is the 
linear span generated by the functions 

Aqx Aia; X n x 
e , e , — , e 

In the case when some Xj occurs rrij times in A n = (Ao, A n ) a basis of the 
space £?(a ,...,a„) is given by the linearly independent functions 

x s e X]X for s — 0,1, mj — 1. 

In the case that A = ... = A n the space -E(a ,...,a„) is just the space of all 
polynomials of degree < n, and we shall refer to this as the polynomial 
case, and we shall denote the fundamental function by $ po i,n ( x )- Obviously 
$poi,n ( x ) is of a very simple form, namely 

%oX,n(x) = -X n , (12) 

so the Taylor expansion is evident from (fl2|) . Moreover one has a very elegant 
recursion, namely 

$ P oi,n+i {x) = • x ■ $ po i jn (x) . 
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Generally, recursion formulas for the fundamental function <3>A n are not known 
or maybe non-existing, except the case that the eigenvalues are equidistant, 
cf. [21]. We emphasize that the important (and easy to prove) formula 



d_ 

dx 



An+l $(A ,..,A n+1 ) (x) = $(A ,..,A„) (x) (13) 



is not a recursion formula from which we may compute 3>(A ,...,A n+ i) (x) from 
$(a ,...,a„) (a;) • 

Later we need the Taylor expansion of the fundamental function $A n 
which is probably a folklore result and not difficult to prove; as definition of 
the fundamental function we take formula (j4]). 

Proposition 3 The function $a„ with A n = (Ao,--.,A n ) satisfies A % (0) 

/or = 0, n — 1, and for k >n the formula 

oo 

$ (aLa„) (0) = E A o°- A ". 

so+...+s n +n=k 

holds. In particular, $^ Aj (0) = 1 and $("^. 1) jAn) (0) = A + ... + A n , and 

oo 

$ (AotL)(°)= E A 0°-^- (14) 
s +...+s„=2 

Proof. Write for z e C with Id > |A,-| 



_L_ = 1.^^.(1) 

z — An z i- hi ^ 3 \z) 



s+1 



Thus we have 

oo oo 



*(Ao,.,A n ) (*) = E - E i X 0--- X n z s 0+ . e Z + n + l dz - 

s =0 s n =0 JI > 

By differentiating one obtains 

OO OO ^ /i ^» 

$ (A ,...,A„) ^ = E E L X 0--- X ^ z s +...+s n +n+l dZ - 

s =0 s n =0 JFr 

For x = the integral is easy to evaluate and the result is proven. 
Let us specialize the last result: 
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Proposition 4 In the case of (Xq, A n ) = A 2s+1 (—1, 1) the following holds 
(0) = 1 and (0) = and $gj 3) (0) = s + 1 

Proof. The equation (0) = 1 is clear; further (0) = A + 

... + A„ = 0. Next we use formula (jT4"|) : split up the integral according to the 
cases s n = 0, 1, 2 : then 

n 

< 2 L, (o) = (o) +E = ^ (o) + a «*!":. 1, ,a„, (0 ■ 

Using this formula for ( £ > [ A( J hl ' ) An-:L ) (0) (instead of $|™^ 2 ' ) ^ (0)) one obtains 

< 2 '.« (0) = *S (0) + A^,«W A „_ d (0) + A.<^ (0) . 

Applied to A n = A 2s+ i (-1,1) and 2s + 1 = n we see that $J^ 1 . ) )An) (0) = 
and $1"^ (0) = A n _i, so we have 

*K',a„, (0) = *S ln _ a) (0) + l. 

By iterating one has $["^. 2) )An) (0) = ^\ n x ~^x n , 2j) (0) + j. Since n = 2s + 1 is 
odd we can put j = s and obtain 

$£+ + i 3) (0) = 8 + (0) = s + 1 

since $(a ) ( x ) — eA ° x an d Ao = 1. ■ 

We finish the section by recalling two standard facts: 

Proposition 5 The function $a„ is real-valued if (Ao,...,A ra ) is equivalent 
to (A , A n ) . 

Proposition 6 If Xq, A n are real then <£>A n (x) > for all x > 0. 

Proof. Since A , X n are real the space Erx ,-,X n ) * s a Chebyshev space 
over K, so $a„ has at most n zeros in R. Since $a„ has exactly n zeros in 0, 
it has no other zeros. By the norming condition $^' ) (0) = 1 it follows that 
^A n (x) > for small x > 0, hence <3>A n (x) > for all x > 0. ■ 

The following result is a simple consequence of the definition of $a„ 
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Proposition 7 If A n = (X ,...,X n ) and c + A n : = (c + A , c + A n ) /or 

some c G C i7ien 

$ C +A„ (x) = e cx $ An (x) . 
If cA n = (cAo, cA n ) /or t/ien 

^cAn («) = \$A n (cx) 
c n 

3 The fundamental function for A2 S +i (A, /i) 

Let T r (£) = re lt for i G [0, 2ir] and fixed r > 0. Let A and \x be two real 
numbers and assume that r > be so large that A and \i are contained in 
the open ball of radius r and center 0. In the following we want to compute 
and analyze the fundamental function 

1 f e xz 

tj - M< * 1= 55l r ( *-A)"'(*-„r' fa - 

By Proposition H it is sufficient to consider the case A = 1 and /x = — 1, so 
we define 

2s+1K ' 2mJ Tr (z-l) s+1 (z + l) s+1 

This integral can be evaluated by the residue theorem, giving two summands 
according to the poles —1 and 1. By a simple substitution argument we see 
that 

<&2s+l (x) = / rrrdz H / Trrdz 

2s+1 1 ; 2m J Tl z s+i ( z + 2 ) s+1 2m J Tl z *+i ( z _ 2 ) s+1 

where T 1 (t) = e lt for t G [0, 2n] . We substitute in the second integral the 
variable z by — z and obtain 



$2s+i (%) = / rrdz / r-rdz. 

2s+1 1 J 2m J Vl z s+i ( z + 2 ) s+1 2m J Tl z *+l ( z + 2 ) s+1 

Now for an integer a we define the polynomials 



1 

2m J Tl 2r«+i f z + 2) 



^ Q (*) = — / ^ ; , ^. + i +a <** (15) 
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and we see that 



<W (x) = e*P°(x)-e-*P»(-x). 



(16) 



Using equality cosh a; = \ (e x + e x ) and sinh (x) = \ (e x — e x ) it is easy to 
see that 

$ 2s +i (x) = cosh (x) [P° (x) - P s ° (-x)] + sinh (x) [P s ° (x) + P° (-x)] . 
For example, it is easy to see that 

$1 (x) = sinhx and $ 3 (x) = - (x cosh x — sinhx) . 

2 



A similar consideration shows that for the function $ 2s (x) defined as 

1 f e xz 
®2s (x) : = / ti dz 

we have 



(x) 



2th 7 ri (* + 2f 2vn 7 Fl (z - 2) 



rfiz + - — : 



^s+l 



dz 



dz 



2m J Tl z^ 1 (z + 2) s ~~ ' 2tii J Tl z s ( z + 2 ) s+1 
= e^PT 1 (x) + e^Pl, (-x) . 

Here for the case s = we use the convention P° 1 (x) := 0. The following is 
straightforward: 

Lemma 8 The functions P s ° (x) are polynomials in the variable x of degree 
s given by the formula 

Proof. From residue theory it is known that 



dz' 



[e xz (z + 2)-*- 1 ] 



z=0 
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The rule of Leibniz gives 



Note that (z + 2) _s_1 = (s - 1) ... (s - k) (z + 2)~'~ 1 ~ k , and clearly 



£ [<*- (z + 2)-'] = e» £ (') (-1)* li±M!,.- ( 2 + 2)— 



fc=0 

Reverse now the summation, and we arrive at 
£ [e- <* + 2)-'] = e" £ (f) (-!)-' (* + 2p-<-> 

fc=0 ^ 7 

Take now z — 0. ■ 

The following is a short list of the first four polynomials for a = : 

(*) = ^ 
P?(x) = \(x-1) 
P°(x) = ^{2x 2 -6x + 6) 

P 3 ° ( x ) = JL ^_20 + 20a; - 8ar 

1 / 20 2 

P° ( x ) = — 70 - 70x + 30x 2 x 3 + -x- 

4 v ; 2 9 V 3 3 

The polynomial P 2 ° (x) is strictly positive on the real line, hence the polyno- 
mials P s ° (x) are not orthogonal polynomials with respect to any measure on 
the real line. However the following is true: 

Theorem 9 The polynomials P s ° satisfy the following recurrence relation: 
As (s + 1) P s ° +1 (x) = x 2 Pl, (x) - 2s (2s + 1) P s ° (x) . (17) 



2 ' -x 3 
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This can be derived by a direct but somewhat tedious calculation. For 
the fundamental function $2s+i — e x P s ° (x) — e~ x P® (— x) we obtain by a 
straightforward calculation the recurrence relation ([TBI below. Since we shall 
derive this recurrence relation from Theorem [28] in Section [5] by a different 
method, we omit the proof of ([1] 



Corollary 10 The fundamental function $2s+i (x) satisfies the recursion 

$2s+3 (x) = X 2 t- r$ 2 s-l (x) ~ ^ + ^ &2s+l (x) , (18) 

2s+i v ) 4s (s + 1) v ; 2s + 2 + v ' y ' 

and the following estimate holds for all x > : 

o< $2s+l( ; T | < / - . (19) 

" $ 2s _! (x) 2s (2s + 1) K J 

Proof. Since $2s+3 (%) > the equation ([TBI implies that 

^^$ 2s+ l (X) < X 2 / $2,-1 (X) 

2s + 2 4s (s + 1) 

from which (fT9|) is immediate. ■ 

Let us recall that after formula ([5]) we defined $2s (%) as the fundamental 
function for the vector with s + 1 many 1 and s many —1. By $2s- (x) we 
denoted the fundamental function with respect to the vector with s + 1 many 
—1 and s many 1. We shall denote sometimes $2s (x) also by $2s+ (x) in order 
to facilitate some formulas and to underline the difference to $2s- (x) ■ The 
following simple identity 

1 1 2 



(z-l) s+1 (z + l) s (z-l) s (z + l) s+1 (z- l) s+1 (z + l) s+1 
implies the formula 

$ 2s (x) - $ 2s - (x) = 2$ 2s+ i (x) . 

This formula can also be derived by summing up the following two identities 
in Theorem [11] which we shall derive from Theorem [271 

Theorem 11 The following two recursions hold: 

$ 2s +i (a;) = $ 2s (x) - —x ■ $ 2s -i (x) , 

$2s+l (X) = -$ 2s - (X) + — X • $2,-1 (x) . 
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Proof. We derive the result from Theorem [23 We choose k = 2s — 1 in 
the relation 

^±P2s+l,2s+l 0) = X ■ p 2s -l,2s-l ~ 2S • p 2s ±,2s (x) 

or which is the same (up to notation) 

A±$ 2s+1 (x) = x ■ $ 2s -i - 2s ■ $ 2s± (x) . 

By Proposition [3] 

^_ 1 (0) = A + .... + A 2a _ 1 = 
$g s ± +1) = A + ... + A 2s± = ±l. 

So A ± = {2s + 1) $£1 (0) - 2 s $g s ± +1) (0) = ±1 (-2.) . ■ 
Corollary 12 T7ie following limit exists 

hm — -> 1. 

Proof. By Theorem [TT1 we have 

$2s (X) _ ®2s+l jx) + Ys X ■ $ 2,-l (g) 
$ 2s _ (x) ^X • $ 2s _! (x) - $ 2s+ l (x) ' 

Let us define y s := $ 2s+ i (x) /$ 2s _i (x) , then 

$2s (a) _ g • Va + |g 
$2 S - (g) \x-s-y s ' 



(20) 



From Theorem [TU1 we see that s ■ y s converges to 0, so (1201) converges to 1. ■ 
We mention that one can derive also recursion formula for the derivatives, 
e.g. the following identity holds: 

Theorem 13 The derivatives J^-P s ° of the polynomials P s ° (x) can be com- 
puted by 
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4 Generating functions 



The Lagrange inversion formula, see e.g. [I], is another way to investigate 
the polynomials defined in (TT5T) 



P° ( x ) = / jAjrdz, (21) 

n 1 ^ 2vu J Vl z n+l ( z + 2 ) n+1+a 1 ; 

where a is a fixed integer. Since this powerful method is somewhat technical, 
let us recall the basic facts. In our case we put tp (z) = l/(z + 2) and 
/' (z) := e xz I (z + 2) a . The fundamental idea of Lagrange inversion is based 
on the observation that 

can be seen as the n-th Taylor coefficient of a holomorphic function which 
will be constructed from /' and (p. More generally, we may assume that <p is 
a holomorphic function in a neighborhood of such that ip (0) ^ 0, and / 
is holomorphic in a neighbhorhood of 0, and we define a n by the expression 
( |22|) . Consider the function 

y(z) :-- 



ip(z) 

which is holomorphic in a neighborhood of (since (f (0) ^ 0) with y (0) = 0. 
Since obviously y(z)ip(z) = z we obtain ip(z)y'(z) +y(z)ip'(z) = 1, so 
y' (0) 7^ 0. Hence y is injective in a neighborhood; let y~ x be the inverse map. 
Since y~ 1 oy (z) = z one has -j- z y~ l {y (z)) ■ -^y (z) = 1, and using this formula 
one arrives at 

(| (/o,-)) M = /V to))- (s^'W^'W'J^)- 

ior y = y (z) . Thus we obtain 



n+1 n ,/'-v^ n + 1 „,f^\ n+1 



y [z) y (z) 



(23) 



Let now 7 (t) = re lt for r > sufficient small, and put z — 7 (i) . Note that 
r (t) = y(j(t)) is a path surrounding zero. Insert z = 7 (i) in f[2"3"j) and 
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multiply it with 7' (t) /2m, so we obtain 

a n = — [ !^L\ v (z)] n+l dz 
2m J z n+1 Jl 



1 F*t y U°y- l ){y{i{t)))-i z y{i(t)) 

' 7 (t) at 



2 «io ?/(7 0O) 

1 fi y Uoy- l ){y) 

-dy. 



n+l 



2m J T y n+1 

Let us define F (y) = ^ (/ o y^ 1 ) (y) , then by the Cauchy formula for the 
Taylor coefficients one obtains 

Since the Taylor series of F (y) is equal to F (y) we obtain 

£ a n y n = F(y) = -f(fo y~ l ) (y) = f (y' 1 (y)) (y) . (24) 

n=0 ay ay 

So this means that for the computation of the unknown sum Y^=o a ny n we 
only have to compute y~ l (y), and the derivative (y) and then we have 
to compute the right hand side of (1241) . As an application we prove: 

Theorem 14 Let (x) be defined by the expression l[21\) . Then for all y 
with \y\ < 1 



Fpa (x,y) :=J2P:(x)y n = e 



-x c 



n=0 



2V^TT(V^TT+1) 



Proof. Let (p (z) — [z + 2) 1 and y (z) := z/cp (z), so y (z) = z (z + 2) 
We define the inverse function y~ l of y (z) by 



y 



^{y) := v^ + T-1. 



Then y 1 (y (z)) = z (z + 2) + 1 — 1 = z and note that 

dy- 1 1 

(y) 



dy 2v^TT' 
17 



Define /' (z) = e xz / (z + 2) a . Then the Lagrange inversion formula tells us 
that 

V P° (x) V n = f ( \/yTl - l) t-== ra 



n=0 

which is exactly our claim 



Remark 15 It is easy to derive a second order differential equation for the 
generating function y 1 — > Fpo (x, y) , and this can be used to give a proof for 
the recursion /TT7D. 



Theorem 16 The following identity 

Fodd 0, y) = ^2 $2s+1 ( x ) ' yS = ^y~+J sinh ( x ^ y + 1 ) 

holds for all y E C 

Proof. Since $2s+i (%) = e x P^ (x) — e~ x P® (—x) we obtain 

00 00 

Fodd (x, y) = e x J2 p s (*) y s - e~ x p s (-*) v*- 

Using Theorem [TH a short computation shows that this is equal to 

F odd (x,y) = -rL=\ ( e *Vv+i_ e -xVv+i\ = _J_ sinh x ^^7. 

v y + 1 2 v y v 2/ + 1 

Since 

00 

£ P n ° (*) y» = exp (x - l)) — _ 



n=0 

we obtain 



Fodd (s,y) = (xJV+i) ^= - exp (-zyfi+l) ^== 
= - ? =L=sinh(x v ^TT). 

Since sinh z contains only odd powers in the Taylor expansion it follows that 
y 1 — ► Fodd {x,y) is entire. ■ 
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Let Ao, Ax, .... be a sequence of real numbers. In the following we consider 
the generating function 



(x, y) := J2 $ (Ao,..,a„) (x) ■ y n 



n=0 



for the case that \2j = 1 and A2J+1 = — 1 for all j e No- With the notation 
from the introduction we have 

$(A f) ,...,A 2s ) = $ 2s +i and ^(Ao,...^.!) = ®2s 
and now we consider 

00 00 
Fa (x, y) = J2 (x) • y 2s+1 + E ( x ) ' V*- ( 25 ) 

s=0 s=0 

It is easy to see that + l) $2s+i (x) = $2s (x) (cf. the general formula 
ffl3]) ). so we have 

E (*) • = + x ) E (*) • y 2s = (l + 1 ) F ^ (*> y 2 ) ■ 

s=0 ^ ' s=0 ^ ' 

Thus we have 

(x, y) = (l + y) F odd (x, y 2 ) + ^^odd (ac, ?/ 2 ) 
and the following is proved: 

Theorem 17 W^i/i i/ie above notations, the generating function F\ (x,y) in 
( ESP zs egna/ to 



F A (x, = 7=== sinh fi^/t/ 2 + l") + cosh (a; a/ y 2 + l") 



5 Construction of Bernstein bases 

In this Section we return to the general theory of exponential polynomials 
where the eigenvalues Ao, A n may be complex numbers. We shall character- 
ize the spaces £Va ,...,A«) of exponential polynomials which admit a Bernstein 
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basis. Let us emphasize that the existence results already follow from those 
in [5], [12], [25] in the more general context of Chebyshev spaces. We 
follow here our exposition in [2] which is based on a recursive definition of the 
Bernstein basis, and it seems that this approach is be different from those in 
the above cited literature. Further references on properties of Bernstein-like 
bases are [TO] and [29|. 

Let us recall some terminology and notations: The k-th derivative of a 
function / is denoted by f^ k \ A function / G C n (J, C) has a zero of order 
k or of multiplicity k at a point a G I if / (a) — ... = f^ 1 ^ (a) = and 
f^ (a) ^ 0. We shall repeatedly use the fact that 

k\ ■ lim -^L = /« (a) . (26) 
x '^° (x — a) 

for any function / G C^(7) with / (a) = ... = (a) = 0. 

Let us recall the definition of a Bernstein basis: 

Definition 18 A system of functions p n ^, k = 0, ...,n in the space Efx ,...,x n ) 
is called Bernstein-like basis for -E(a ,...,a„) ond a ^ b if and only if each 
function p n ± has a zero of exact order k at a and a zero of exact order n — k 
at b. 

In the following we want to characterize those spaces Er\o,-,\ n ) which 
admit a Bernstein basis. For this, we need the following simple 

Lemma 19 Suppose that fo, /„ G -E(A ,...,A n ) have the property that fk has 
a zero of order k at a for k = 0, n. Then f Q , ...,/„ is a basis for any interval 
[a, b] . 

Proof. It suffices to show that f ,...,f n are linearly independent since 
^(Ao,...,An) nas dimension n + 1. Suppose that there exist complex numbers 
Co, c n such that 

Co/o (x) + .... + C n f n (x) = 

for all x G [a, b] . Since fo (a) 7^ and fj (a) = for all j > 1 we obtain 
Co = 0. Next we take the derivative and obtain the equation 

c 1 f[(x) + .... + c n f' n (x)=0. 

We insert x = a and obtain that c\ — since f[ (a) 7^ and /j (a) = for all 
j = 2, ...,n. Now one proceeds inductively. ■ 
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It follows from the Lemma [TH] that a Bernstein basis p Ut k,k = 0, ...,n 
is necessarily a basis for Er\ 0t „. t \ n y Next we want to show that the basis 
functions are unique up to a factor (provided there exists such a basis). 

Proposition 20 Suppose that there exists a Bernstein basis p n ,k, k = 0, ...,n 
for the space i?(A ,...,A„) an d a,b. If f G -E(a ,...,a„) has a zero of order at least 
k in a and of order at least n — k at b then there exists a complex number 
c such that 

f C • Pn,ko ■ 

Proof. Since p n ,k, k = 0, ...,n is a basis we can find complex numbers 
Co, c n such that 

/ = CoPnfl + ■■■ + C n p n , n . 

We know that / (a) = ... = f^ ^ (a) = since / has an order of exact order 
k in a. So we see that = / (a) = c p n fl (a) . Since p n $ (a) ^ we obtain 
c = 0. We proceed inductively and obtain that 

f CkoPn,ko CnPn,n- 

Now we use the zeros at b : Inserting x = b yields 

= / (6) = c n p ntn (b) . 

Since p n , n (6) ^ (here we use the exact order at the point b) it follows that 
c n = 0. Proceeding inductively one obtains / = Ck p n ,k - ■ 
The proof actually shows the following: 

Corollary 21 Suppose that there exists a Bernstein basis p n ,k, k = 0,...,n 
for the space -E(a ,...,a„) and a, b. If f 6 E(\ ,...,x n ) has more than n zeros in 
{a, b} then f = 0. 

Let us recall the following definition: 

Definition 22 The space E(x ,...,\ n ) is a Chebyshev system with respect to 
the set A C C if each function f 6 ^(Ao,...^™) that has more than n zeros in 
A (including multiplicities) is 0. 

So we have seen that a necessary condition for the existence of a Bernstein 
basis on the interval [a, b] is the property that -E(A ,...,A n ) is a Chebyshev system 
for the set {a, b} . In the next Proposition we shall show that this property 
is actually equivalent to the existence. 
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Proposition 23 The space E^x ,...,x n ) is an extended Chebyshev system with 
respect to a, b if and only if there exists a (unique) Bernstein-like basis p n ± '■ 
k — 0, n, for a, b satisfying the condition 

k\ lim f^=P { n k i(a) = l. (27) 

The Bernstein-like basis functions p n ^ (x) are recursively defined by equations 
KM), jm fffl and M) below. 



Proof. The necessity was already proved. Assume now that E(x 0r ..,\ n ) is 
an extended Chebyshev system for a, b. It is convenient to use the following 
notation: 

q k (x) := p n ,n- k (x) , (28) 
so qk has a zero of order n — k at a and a zero of order k at b. Define first 

go (x) := $ An (x-a), (29) 

which clearly has a zero of order n at a and of order at least at b. Since 
go can not have more than n zeros on a, b by our assumption we infer that 
g has a zero of order at b. Define gi := g^ — a go for a = q^ (b) /g (b) 
which has a zero of order at least n — 1 at a and a zero of order at least 1 at 
b. Again our assumption implies that gi has a zero of order at n — 1 at a and 
a zero of order 1 at b. For k > 2 we define q k recursively by 

qk ■= g£_i - («fc-i - a fc _ 2 ) • gfc-i - (3 k qk-2 (30) 

with coefficients a k -i — a. k -2 and (3 k to be determined. By construction we 
know already that gfc„i and q k -2 respectively have a zero of order k — 1 and 
k — 2 at b, and a zero of order n — k + 1 and n — k + 2 at a. So it is clear that 
q k has a zero of order at least k — 2 at b, and a zero of order at least n — k at 
a. The coefficients — and f3 k are chosen in such a way that g^ will 
have a zero of order at least k in 6. This is achieved by defining 

P* '= (fc-2) aIld := ( 31 ) 

Qk-2 ( b ) Qk-1 ( b ) 

Again our assumption implies that q k has a zero of order k at b and a zero 
of order ra — k at a. The condition f|2T|) is easily checked using (|2"fi]l . and (I3H 
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together with induction. The uniqueness property follows from the above 
remarks. ■ 

In the rest of the paper we shall call the Bernstein-like basis provided by 
Proposition [231 the Bernstein basis with respect to a,b. 

Next we shall give a construction of the Bernstein basis p n ,n-k 0*0 , k — 
0, ...,n, will is similar to constructions known from the theory of Chebyshev 
spaces. 

Theorem 24 Let (A ,...,A n ) G C n+1 and define for each k = 0,...,n the 

(k + 1) x (k + 1) matrix A n<k (x) by 



A n k (x) 



**(x) ... *<P(x) \ 



*<*>(*) ... < fe) (*)/ 



(32) 



Then the matrices A n ^ (b — a) are invertible for k — 0, ...,n if and only if 
E(\ 0t ...,\ n ) is an extended Chebyshev system with respect to {a, b} . 

Proof. Assume that the matrices A n ^ {b — a) are invertible for k = 
0, n. It suffices to show that there exists a Bernstein basis with respect to 
a, b. For a polynomial r?y n ^ (z) = r k z h + r k _ 1 z k ~ 1 + ... + r with coefficients 
r , r k G C and r k = 1 let us define 

U (x) :=±[ rAn '"[ Z) f~ a \ S dz. (33) 
2m J (z-Xo)... {z- A n ) 

It is easy to see that f n ,k is in £7(a ,...,a„) since 

k 

fn,k (x) =5^r i $S2(x-o) 

3=0 

From this representation it follows that has of order at n—k at a, and that 
j(n-k) ^ _ rjfc $w ^ x _ a ^ = i ^y e wan ^ choose the coefficients ro, r k -i 
such that ^ 

■^ifn,k) -I" for all Z I) /,• I 
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and that \J^f n ,k) (b) 7^ 0. Writing down these equations for / = 0, k 
in matrix form shows that for c = b — a 



(fe-1) 

A n 



(fc-1) 



(2fc-2) 



r 



\ 



(c) / V r *-i / 



\ ^A„ 



S(c) 



/S ( & ) = o- 



Since the matrix A n ^_i (b — a) is invertible we can find clearly r , r^-i G C 
which solves the equation. Hence f n ^ has a zero of order at n — k at a and 
a zero of order at least k in b. Suppose now that / (b) 
Then these equations say that 

^ ^ f 



(2k) 



)l 


^r \ 







r fc = 0. 



Since the matrix A n ^ (b — a) is invertible it follows that r = ... 
This is a contradiction to the choice = 1. 

Now assume that £?(a ,...,A») is an extended Chebyshev system with respect 
to {a, 6}. Suppose that the matrix A n ^ (b — a) is not invertible for some k G 
{0, n} . Then we can find s = (sq, Sjt) 7^ such that A n ^ {b — a) s = 0. 
Define / n> fc (x) = Y^j=o s j^aI ( x ~ a ) ■ Then clearly f n ^ has a zero of order 
at least n — k at a. Further the equation A n ^ (b — a) s = implies that 
fnl (b) = for j = 0, k. So / n> fe has a zero of order n + 1 in the set 
{a, b} , and our assumption implies that f n ^ = 0. By Lemma [191 the system 

$A n ,---,$A is a basis of -E7(a ,...,a„)- It follows that so = ... = 
contradiction. ■ 

The proof of the last theorem actually shows: 

Theorem 25 Assume that the matrices A n ^ (b — a) are invertible for k 



0, a 



1, n + 1, and let k = 0, 

a, 6. T/ien /or eac/i k = 0, ...,n there exists a polynomial r^ njk (z 
and leading coefficient 1 such that 



n be the Bernstein basis with respect to 

of degree k 



Pn,n—k (•^') 



1 

2tH 



riv n ,k{z)e {x a)z 
{z - A ) ... (z - \ n ) 



dz 



for all x G 



24 



For a given vector A n = (A , A n ) the matrices A nk (z) ,k — 0, n are 
defined in ( 152"j) and we set 

n 

Za„ ■= \J{zeC:detA n , k {z) = 0} . 

fc=0 

Note that Z^ n is a discrete subset of C since z i — > det A„ )fe (z) is obviously 
an entire function. It follows that for given a e M the space E(x ,...,x n ) is 
an extended Chebyshev system for {a, b} for all & E R except a countable 
discrete subset of R. We emphasize that this does not imply that E(x ,...,\ n ) is 
an extended Chebyshev system for the interval [a,b] . In order to have nice 
properties of the basic function p n fc one need the assumption that £7(a ,...,a„) 
is an extended Chebyshev system for the interval [a, b] . Indeed, the following 
result is well known (at least in the polynomial case): 

Theorem 26 Suppose that E(x ,...,\ n ) is an extended Chebyshev system for 
the interval [a,b] with Bernstein basis functions p n ^,k = 0, ...,n. Assume 
that E/x ,...,x n ) is closed under complex conjugation. Then the basis functions 
Pn,k{x) are strictly positive on (a,b) for each k = 0,...,n. For each k = 
1, n — 1 there exists a unique a <t k < b such that 

Pn^k is increasing on [a, tk] 
Pn t k is decreasing on [tk, b] . 

So p nt k has exactly one relative maximum. The function p n ^ n is either in- 
creasing or there exists to G (a, b) such that p n ,n is increasing on [a, t ] and 
decreasing on [to, b] . 

Proof. 1. Consider (x) = p U)k (x) -p n> k (x) for k = 0, n. Then f Utk 
has a zero of order at least k + 1 in a and a zero of order at least n — k in b. 
Moerover / e E(x ,...,x n ) since -E(A ,...,A n ) is closed under complex conjugation. 
As -E(A ,...,A n ) is an extended Chebyshev system for {a, b} we infer f n k = 0, 
hence p n)k is real-valued. 

2. Since £7(A ,...,A n ) is an extended Chebyshev system for the interval [a, b] 
it follows that p n ^ has no zeros in the open interval (a, b) . By the norming 
condition (a) = 1 it follows that p n ^ is positive on (a, b) . Since the 
derivative p' nk is again in -E(A ,...,A n ) it follows that p' nk has at most n zeros 
on [a, b] . 
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3. We assume that 1 < k < n — 1. Then p' nk has a zero order k — 1 > at 
a and a zero of order n — k — 1 at b. Hence p' n k has at most 2 zeros in the open 
interval (a, 6) . First assume that p' n k has two different zeros to and t\. Then 
they must be simple. It follows that p' n k (t) > for t e [a, t ] (recall that p n ^ 
is positive, so it must increase at first). For t E [to,ti] we have p' nk (t) < 
since the zeros are simple. Finally we have p' nk (t) > for t e [ti, fe] . So 
p n ^k is increasing on [ti, b] , which implies p nj fe (t\) < p n ^ (b) = (recall that 
k < n — 1). So p ni fe has an additional zero in t\ which is a contradiction. 
Now assume that p' n k has a double zero at t - Since it has no further zeros 
P'n k(t) > f° r a ^ * e ( a ' ' so is monotone increasing, which gives a 
contradiction to the fact that p n ^ (b) = 0. 

In the next case we assume that p' n k has no zero in (a, b) . Then p Ut k is 
strictly increasing, so p n ^ ih) > 0, which gives a contradiction. So we see 
that p' n k has exactly one zero in (a, b) . 

5. Since n has a zero order n — 1 at a, it has at most one zero in (a, 6) . 
If p' n n has a zero then p n ^ n is increasing on [a, t ] and decreasing on [to, b] . If 
p' n n has no zero then p n>n is increasing. ■ 

We mention that -E(a ,...,a„) is an extended Chebyshev system over inter- 
vals [a, b] whose length b — a is sufficiently small. Moreover, for real eigen- 
values Ao,...,A n the space E^ ,...,\ n ) is an extended Chebyshev system over 
any interval [a, b] . 

Simple experiments show that the assumption of a Chebyshev system 
over the interval [a, b] are essential. In the case of the six eigenvalues 

±7i, ±i (7 - tt) , ±i. 

we have found basis function p n>k for the interval [0, 3] or [0, 3.14] with several 
relative maxima. Even it may be possible that p n ,k are non-negative (but they 
are always real- valued) . The fundamental function </? 6 (x) has its first zero at 
3.2. 

6 Recurrence relations of the Bernstein basis 

Let us recall that A 2s+ i := (s + 1 times A, s + 1 times p) and let us introduce 
the following notations 

A 2s+ := (s + 1 times A, s times p) , 
A 2s _ := (s times A, s + 1 times p) . 
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In this section we want to give recursive formula for the Bernstein basis 



P2s+l,k '■ — Pk 2 s+i.k and P2s+,k '■ — PA 2s +,k- 

Theorem 27 For each k = 0, ...,2s — 1 the following recursion formula holds 

A±p 2 s+l,k+2 (x)=X- p 2s -l,k ~ {k + 1) • P2s±,k+1 (X) 

where the constant A± is given through 

A ± = {k + 2) pf s % (0) - (k + 1) pj&J+i (0) • 

Proof. Let us recall that px„,k is the unique element in E\ n which has k 
zeros in and n — k zeros in 1 with the normalization 

ajSS*^ =*£!«» = i- Pi) 

Consider f±(x) = a ■ x ■ P2s-2,k + & • P2s±,fc+i (x) for coefficients a and b which 
we want to define later. Clearly f± is in E\ 2g+1 . Note that P2s-i,k has a zero 
of order 2s — 1 — in l,and that P2s±,fc+i (#) has a zero of order 2s — (k + 1) 
in 1. Hence f± has a zero of order at least 2s — 1 — k in 1 (just as p2s+i,fc+2)- 
For similar reasons we see that f± has a zero of order at least k + 1 in 0. We 
choose now the constants a and 6 in such a way that f± has a zero of order 
k + 2 in and that / satisfies the normalization in By the uniqueness 
we infer that /± = P2s+i,jfc+2- Clearly /± has a zero of order k + 2 in if 

fl ' (^3TT ( x ' P*-V*]j ^ + & • PaSS+i (°) = °- 
Recall that p^i'fc+i (0) = 1 anc ^ 

/ \ i k , pL-i,*; (o) fc+i . pL-i,fc (o) k+i , 



It follows that 

\ J k+l r i . . . i 

= k + 1. 



fc! 



x=0 
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Thus b = —a (k + 1) . The normalization condition gives the second equation 
1 = vtXtk+2 (0) = a • (jl^ [x ■ P2,-i*]j ~a-(k + l) p£± 2 fc+i (°) • 



Again we see that 



(/ " ' " [x ■ p 2s -i,k (x)} x=Q = (k + 2) p%% (0) • 



s * = £rf + x ) ^-?, +2 (o) - (* + 3) vl% (o)) 



So we obtain 

i = a ((* + 2) P ^;; fe (o) - (* + (o) 

The proof is complete. ■ 

Theorem 28 For each k = 0, 2s — 3 holds the following recursion formula 
A k ■ P2s+l,k+4 = x 2 ■ p 2s -3,k -(k + l)(k + 2) p 2s -i,k+2 + B k -x- p 2s -l,k+2 (35) 
where the constants A k and B k are defined by 

k + 2 
k + 3 

A k = {k + 3) (k + 4) pl% (0) — (k + l)(k + 2) p£%. a (0) 
+ (fc + 4)pg%. a (0).S*. 

Proof. Consider / (x) = a ■ x 2 ■ p 2s -3,k + b ■ p 2s -i,k+2 + c-x- p 2s -i,fc+2 for 
coefficients a, b and c which we want to define later. Clearly / is in E\ 2s+l . 
Note that the functions P2s-3,k,P2s-i,k+2 and p2s-i,k+2 have a zero of order at 
least 2s — 3 — in 1. Hence / has a zero of order at least 2s — 3 — k in 1 (just 
as P2s+i,k+i)- Clearly / has a zero of order at least k + 2 in 0. We choose now 
the constants a, and c in such a way that / has a zero of order k + 4 in 
and that / satisfies the normalization in fl34l) . By the uniqueness we infer 
that / = P2s+i,&+4- Clearly / has a zero of order k + 4 in if and only if 

f(k+2) (Q) = /( *+3) (Q) = Q 

Recall that P 2 s-i\+2 W = 1- Since x -p2s-i,fc+2 has a zero of order k + 3 in 
the equation /( fc + 2 ) (0) = is equivalent to 



/ d k+2 \ 

a - [x 2 -P2s-3,k] ) x=o + b = - 
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-3,fe 



,k+2 



Consider the Taylor expansions of p 2s - 

_ 1 „fc+2 , ^2s-l,fc+2 V U ^ fc+3 , 

' 2s - 1 ' fc+2 " (fc+~2)f + (fc + 3)! * + 



(x) and p 2s -i,i 

P2s-3,k K U ) fc+2 , 
^ _i_ OM ••• 



It follows that 

= -a{k + l) (k 
( d k+3 \ 

Clearly 



a 



\ r I "I 

3, fc (x)] =^ -^+ 2 = (* + !)(* + 2). 
/ x =o ux L"'- Ji=o 

Hence 6 = -a (k + 1) (A; + 2) . The equation f (k+3 *> (0) = implies that 

' P2 S -3,fe]J +^ 2 i+?; fc+2 + C^^ j ^[x-p 2 .-l,fc+2]J =( 



and similarly 



= / d* +3 P^3? fc (0) fe+ 3 N 

~~ I cfc fc + 3 (k + l)\ X . 
= (* + 3) (k + 2)p£%(0) 



ith b = — a (k + 1) (k + 2) we obtain the equation 
+ 3) (k + 2) ^ (0) - a (k + 1) (k + 2) p£t + 2 (°) + c(/c + 3) = 0, 



So with 6 

a (A; + 3) (fc + 2)p£J| fc (0) 
which shows that 

c (fc + 3) = a (k + 2) ((fc + 1) p£-t + 2 (0) - (* + 3) p^Sfc («)) 
so c = Bj,. The normalization condition gives the third equati< 

1 = p£^ fc+4 (0) = a • [x 2 • p 2s _3, fc (x)] ^ ^ 

-!)(* + 2) pl + _% +2 (0) + c ^ [x ■ p 2s -i, k+ 2 (*)] J ^ . 



S,fc 1 

tion 



-a(fc + : 
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Again we see that 
/ d k+A 



and 



\ aX / x=0 



P2s-3,k \ V ) k 



-X 



ydx k + 4 I (fc + 2) 
(* + 3) (* + 4)p£3> (0) 



i=0 



/ d fc + 4 . 



s-l,fc+2 



(*)] 



(k + 3)! 
(* + 4)^^(0). 



/>'+! 



x=0 



So we arrive at the equation 

a- 1 = (k + 3) (* + 4) p^Jj, (0) - (k + 1) (* + 2) ^t!,+ 2 (0) 

(fc+3) 



+ (* + 4)jCii+2(0)S*- 



Corollary 29 Suppose that /x = —A. Then the fundamental function &2s+i (x) 
satisfies the following recursion 



1 2s — 1 

X 2 • $ 2s -3 0) $2,-1 (s) 



4s (s - 1) 



2s 



Proof. Let us take = 2s — 3 in equation (1331) : then p 2 s+i,2s+i = $25+1 
and p 2s -3,2s-3 = *2 S -3 and p 2a -i,2s-i = $2s-i- By Proposition [3] 



(2s-2) 
2s-3 



(0) = (s-2)A + (s-2)/i = 0. 



and similarly $2s-i = 0- Hence for = 2s — 3 the last summand in (1331) is 
zero and the formula amounts to 

A 2s -3 • $ 2s +i = x 2 ■ $ 2s _3 - (2s - 2) (2s - 1) $ 2s _!. 

We have to compute the constant for fc = 2s — 3. Since $2s-i (0) = anc ^ 
$ 2 s-3 2 ' ) (0) = the constant is zero and we have 



A 2s _ 3 = 2s (2s + 1) $gr 3 1} (0) - (2s - 2) (2s - 1) <&£^ (0) 



(2s+l) 
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Proposition H] shows that ^s+i 3 '' (0) = s + 1, so $2s-3 (0) = s — 1 and we 
see that 

A 2s - 3 = 2s (2s + 1) (s - 1) - (2s - 2) (2s - 1) s 
= 2s (s - 1) (2s + 1 - 2s + 1) = 4s (s - 1) . 

Hence 

+ 4s (s - 1) 2s 

■ 

Let us use the notation q 2s +i,i = P2s+i,2s+i-«- Then the recurrence relation 
in this notation means (with I = 2s + 1 — (k + 4), so k = 2s — 3 — /) that 

A 2s -3-rq 2 s+i,i = x 2 -q 2s ~3,i-(2s - 2 - I) (2s - 1 - I) q2s-i,i J rB 2s ^ r x-q 2s -i t u 
so we have a recurrence relation for fixed I. 



7 The derivative of the Bernstein operator 

Let C N ° be the set of all sequences y = (yo,yx,....) = (yi~) with complex 
entries. For y 6 C N ° we shall use the notation 

(y)k '■= Vk 

for describing the fc-th component of the vector y. The finite difference op- 
erator A : C No C No is defined for y = (y Q , y x , ....) = (y k ) G C N ° by 

(^y)k := Vk+i ~ Vk- 

Higher differences are defined inductively by setting A n+1 y = A (A n y) where 
A is defined as the identity operator. The difference operator A is useful in 
the classical theory to describe the derivative of the Bernstein polynomial 



B n f(x) 



k=0 



k /i \n—k 
X (1 — X) 



namely 



n 



k=0 



f\ k - 



n 



n — 1 
k 



x k (1 - x) 



n— 1—k 
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We want to derive an analog for the Bernstein exponential polynomial. 
We recall from [2] the following numbers 

d k /> := hm ^ k {X) 



x ^ b PA„\\j,k (x) 

which have been important for the construction of the Bernstein operators. 
We shall assume that Ao = which facilitates the formulas. We define a 
difference operator A A „ >A . : C N ° -> C N ° for y = (y k ) G C N ° by 

{A An>x .y) k := 4^ • Vk ~ 4'„ A ° • Vk+i 
By A„ \ Xj we denote the vector where we have deleted Xj . Now we can prove 

Theorem 30 Let Xq = and j G {0, 1, ...,n} . For -B(a ....,a„)/ the following 
identity holds: 

d n ~ 1 

Proof. From the definition of the Bernstein operator we immediately see 
that 

d n d 

- ^j) B (\a,-,K)f {x) = ^ (fa ~ A j)P(Ao,..,A„),fc- 

fc=0 



By the next Theorem below we obtain 

d U n—l 

(~f- ~ ^j) B (Xo,...,X n )f ( x ) = a kf (t k ) PA n \\j,k-l + Yl "*/ d AnPA n \^,k 



1 dx 

k=l k=0 



n-1 

^ [a*/ {t k ) d k ^ 3 + a k+1 f (t k +i)] ■ p.\„ a,./, 



k=0 

From the construction of the Bernstein operator in [2] the following formula 
is known: 

Otk+1 — ~ a kO'A n ■ 

The proof is accomplished by identity 

a k+1 f (t fe+ i) + a k f (t k ) 4' Aj = afc -<fi£°f (tk+i) + f (t k ) dS 



From [2] we repeat the following result: 
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Proposition 31 Define for k = 0, n — 1 



d k A X - . > := lim ^ '-'^^^ (36) 
(A0 '-' Aw) *T6 PCAo,..,^.!),* {x) 

Then, for any k = 1, ...,n — 1 , 

_ P ( Ao '-' A «)' fe = P(Ao,...,A„_i),fc-l + d (Ao".,A„)^o,...,A„-i),fe- (37) 

Furthermore, for k = we have 

-j- ~ An j P(Ao,...,An),0 = rf (A '!!..,A„)P(Ao,.-,An-i),0, (38) 



w/iz/e /or k = n, 

d 



dx ~ ^"J P ( A °'-' A ")' n ~~ P(Ao,...,An_i),n-l- (39) 

Proof. We may assume that a = 0. Set 

/fc := _ ^nj P(Ao,...,A„),A> 

and let 1 < k < n— 1. Using the fact that /& has a zero of order fc— 1 at and of 
order n-k-l&tb, it is easy to see that f k = c fc p(A ,...,A„_ 1 ),n + 4p(A ,...,A„_ 1 ),fc 
for some constants Ck and d k - We want to show that c k = 1. Note that 
P(Ao,...,a„),a has a zero of order in 0, so 

r r Id , x P(Ao,-,A„),fe (°) 1 

IS? -P=5" = i*=r^.-^* (X) = (fc-l)! = (* = !)! 

where the second equality follows from (1261) applied to ^ * fe (x) and 

([27]). On the other hand, the equation f k = c k p(\ 0l ...,x n ^ 1 ) l k-i + 4p(A ,...,A n _ 1 ),fc 
shows that 



xj.o x^ 1 xio x^ 1 (A; - 1)! 
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using again Proposition [231 (E7|) . Hence c k — \. Next we divide fk by 

(6 — x) n_fc 1 to obtain 

4P(Ao,...,A»-i),fe (g) = .. /fe(s) 



P(A ,...,A n ),fc (x) - A n p (A0) ... )A „) !ifc (x) 



lim , . 

lim i> P ^-^ k ^ 



C T& (b — x) 



n—k—l 



The case k — follows by noticing that / = c?oP(Ao,...,An-i),0) solving for d Q and 
taking the limit as x | b, while the case k = n is an immediate consequence 
of the fact that P(A ,...,A„),n = $A„- ■ 
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